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Abstract 

Elastic nN scattering and the reaction n + p — > K + Y, + are described simultaneously in a unitary 
coupled-channels approach which respects analyticity. SU(3) flavor symmetry is used to relate 
the t- and u- channel exchanges that drive the meson-baryon interaction in the different channels. 
Angular distributions, polarizations, and spin-rotation parameters are compared with available 
experimental data. The pole structure of the amplitudes is extracted from the analytic continua- 
tion. 
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1. Introduction 

The excitation spectrum of baryons and mesons is expected to reveal important information 
on the mechanism of confinement as well as the intrinsic structure of hadrons. Properties of 
baryon resonances have been obtained by lattice calculations (llSJ3J3l> mostly for the ground 
states but also for some excited states J2, Et]. In quark models T|ff la, 01, a rich spectrum of 
excited states is predicted. Many of these resonances could be identified in elastic nN scattering, 
while at higher energies, usually more states are predicted than seen, a fact commonly referred 
to as the "missing resonance problem" Since resonances not seen in the nN channel might 
predominantly couple to other channels, there are intensive experimental efforts to measure, 
among others, multi-pion or KY final states, where KY = KA or KI.. 

The reaction n + p — > K + l. + provides access to a pure isospin 7 = 3/2 two-body reaction 
channel in meson-nucleon dynamics. Moreover, the weak decay £ + — > pn° allows to determine 
the polarization of the produced E + . In the 1980's, Candlin et al. measured differential cross sec- 
tions and polarizations at the Rutherford Appleton Laboratory for pion beam momenta ranging 
from p K = 1.282 GeV/c to p n = 2.473 GeV/c d and performed an energy-dependent isobar 
analysis for invariant collision energies ranging from the K + 2, + threshold (z = -\fs — 1.68 GeV) 



to V« = 2.35 GeV ll ill . While the quality of the fit is good, unitarity is violated and a sepa- 



ration of the resonant part is difficult due to the oversimplified construction of the partial wave 
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amplitudes. The resonance parameters extracted confirmed four A-resonances found previously 
in the part ial wave analyses of elastic pion-nucleon scattering by Cutkosky et al. H2I1 and Hohler 
et al. MB, the A(1905)F 35 , A(1920)P 3 3, A(1950)F37, and A(2200)G 37 . Other resonances 
deduced from elastic pion-nucleon scattering by Refs. 1 12, T3L [T3| could not be unambiguously 
identified in K + Y. + production ll 111 , notably the A(1910)P3i, which is given a 4-star status by the 
PDG, and the A( 1900)5 31, downsized to a one star rating (nowadays, two star) after the KI. data 
became available. 

In 1988, Candlin et al. IU5I1 measured the spin-rotation parameter B at the CERN-SPS, us- 
ing a polarized frozen spin target and the Rutherford Multiparticle Spectrometer RMS, adding 
independent information to the data base for two pion beam momenta p„ = 1 .69 GeV/c and 
p n = 1.88 GeV/c. The spin-rotation parameter correlates the spin of the target proton and the 
spin of the produced I + and allows to eliminate ambiguities in the partial wave analysis. Dis- 
crepancies between the /^-values predicted from the partial wave analysis of Ref. 1 1 ill and the 
experimental ones were found which suggested the necessity of a new partial wave analysis 115 1. 

A consolidated knowledge of coupled-channels meson-baryon (MB) scattering is required 
when studying meson production in AW collisions, such as measured at COSY/JulichjLU 
There, the MB — » MB transitions enter the proton induced strangeness production IK 
sub-processes in on-shell but also off-shell kinematics. 

Also, a detailed knowledge of the resonance content in KI production is needed in heavy ion 
collisions. The K + has a long mean free path in the nucleus and is believed to provide information 
about the high density and temperature phase of the heavy ion collision 11911 : to clarify the role 
of the A( 1920)^33 in the n + p — > K + I + reaction is thus mandatory 111 911 . 

Various analyses of meson-baryon scattering are available, designed with the goal to extract 
resonance properties from data. They differ, e.g., in the number of channels and their analytic 
properties. Some representative analyses are discussed in the following. 

A coupled reaction channel analysis of nucleon resonances including the KI channel has 



been performed by the GieBen group in the ^-matrix approximation 11201 12111 . Elastic nN scatter- 
ing, as well as the family of kN — > KY, r]N, 10N and other reactions are included in the analysis. 
Photon-induced reactions within the model have been studied in Ref. 02 111 . The non-resonant part 
of the amplitude is treated in a Lagrangia n appro ach and resonances are included up to a total 
spin of J — 3/2. In more recent studies 11221 12311 . the spin 5/2 resonances have been included 
within a Lagrangian-based framework. Unitarity is respected, but the real, dispersive parts of the 
two-body intermediate states are neglected, such that analyticity is lost. 

While in this analysis the imaginary part from phase space is cut off at threshold, in other 
approaches it is analytically continued below threshold, but the dispersive parts are still not in- 
cluded. Such /T-matrix approaches 12-4 12511 analyze jiN, 77^ or ttttN data, or even more reactions 
like the Bonn-Gatchina group 

A very precise analysis of elastic nN scattering is provided by the ^-matrix approach of 

" HE, 



3211 . There are no assumptions made about resonances [ex- 



the GWU/SAID group J29 
cept for the A(1232)] and in this sense the extraction of the excited baryon spectrum is model- 
independent. Also, this partial wave analysis provides the lowest \ 2 of the available analyses 
of elastic ttN scattering 13 ill . This is also the reason, why in this study we use the analysis of 
Ref. IB ill as input rather than directly fitting to nN data, although a direct fit to data should be 
carried out in the future. Interestingly, in the most recent update of the analysis 13 ill , several 
resonances with less than four stars could not be confirmed any more. 

Carnegy-Mellon-Berkeley (CMB) type of models I33ll34il35ll usually include the dispersive 
parts of the resonance propagators but do not provide a microscopical background. 
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Dynamical coupled-channels models take the real, dispersive parts of the intermediate states 
int o account and provide a microscopical description of the background 

H EE M M M 



41l 1421 143L I44L 14511 . Dynamical coupled-channels models are based on effective Lagrangians. 



While jtN scatte ring at low energies is completely understood from chiral perturbation theory 
(see, e.g., Refs. |4fJ 47] or Refs. (37, 48] for unitarized extensions of ;fPT), at higher energies 
model assumptions need to be made. It is realistic to assume that the interaction is driven by the 
exchange of known mesons and baryons. The scattering amplitudes are then obtained as solutions 
of a Lippmann-Schwinger equation which guarantees unitarity. Thus, the driving term of the 
Lippmann Schwinger equation consists of f-channel meson exchange processes and w-channel 
baryon exchanges as well as s-channel processes which may be considered as bare resonances. 

The explicit treatment of the f-channel and M-channel diagrams introduces strong correla- 
tions between the different partial waves and may generate a non-trivial energy and angular 
dependence of the observables. The explicit treatment of this background in terms of exchange 
diagrams also allows to link different reactions such as elastic nN scattering and the reaction 
n + p — * K + Y. + , using SU(3) flavor symmetry. Thus, the treatment of the interaction via meson 
and baryon exchange is expected to lead to a realistic background, with strong restrictions on the 
free parameters. 

In view of this, the strategy to perform baryon spectroscopy is to introduce only a minimum 
number of bare resonance states in order to obtain a good description of the data. This distin- 
guishes the ansatz from some /^-matrix approaches where the absence of a structured background 
may require the introduction of additional resonance states, which improve the x 1 but are in fact 
simulating the background. 

Dynamical coupled-channels app roaches have been so far restricted to the reaction channels 
Nn,Nr],Ncr, An, and Np |39l 40[ 44], and concentrated on differential cross sections, mostly of 
elastic nN scattering. In the present study, we extend the dynamical coupled-channels Jiilich 
model, which has been developed over the years J38i [39 , 4(| 41 ] , to the kaon-hyperon sector by 
adding Lagrangians for the coupling s to the kaon hyperon channels and resonances beyond the 



set considered in Re fs. 08|. l39L 
investigated in Refs. O8[l39Ll40l 



40 , 14111 . We limit our resonance analysis to the energy range 
4lh . i.e. 2 GeV, and concentrate on the isospin 1 = 3/2 sector. 



In Sec. l2.3l the inclusion of the KY channels in addition to the channels Nn,Nr], Ncr, An, and 
Np is discussed. To describe the data in the n + p — > K + I, + reaction, we also need to include higher 
spin resonances up to a total spin of J = 7/2 (cf. Sec. 12.21 ). Results are presented in Sec. [3] For 
the analysis of the resonance content of the resulting amplitude, given by the pole positions and 
residues, one needs the analytic continuation, summarized in Sec. 12.41 The extracted resonance 
properties are listed and commented on in Secs. l4.ll and l4.2l In |Appendix A| ( [Appendix ~B\ , the 
f- and u- (s-)channel processes used in this study are explicitly given. 



2. Formalism 



2.1. Scattering equation 

The coupled-channels scattering equation 1 38 , 39[ 4^, 4^] used in the present formalism ful- 
fills two-body unitarity, as well as some requirements of three-body unitarity following Ref. iBoll . 
Furthermore, it fulfills analyticity and takes into account the dispersive parts of the intermedi- 
ate states as well as the off-shell behavior dictated by the interaction Lagrangians. This integral 
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equation which is solved in the JLS -basis is given by 



{L'S'k'\T'J y \LSk) = (L'S'k'W'^LSk) 

4 J y \L"S"k , 



+ 2 f k" 2 dk"(L'S'k'\V' u J y \L"S"k") - 1^ : - (L"S"k"\T"\LS 

yL»S>> J Q 



where J(L) is the total angular (orbital angular) momentum, S {I) is the total spin (isospin), 
k(k', k") are the incoming (outgoing, intermediate) momenta, and fi, v, y are channel indices. 
The incoming and outgoing momenta can be on- or off-shell. In Eq. dj]), E y is the on-mass shell 

energy in channel y, E = yjm 2 + (k") 2 + ^m 2 B + (k") 2 where m (»!«) is the meson (baryon) 
mass. The second term in Eq. ([1]) on the right-hand side involves also a sum over all intermediate 
possible quantum numbers and channels contained in the model. 

For the channels involving quasi-particles, <tN, pN, and nA, the propagator is slightly more 
complicated 138L l39l pill (cf. also Sec. \2.4\ . The pseudo-potential V iterated in Eq. dp is 
constructed from an effective interaction based on the Lagrangians of Wess and Zumino I52LI53I1 . 
supplemented by additional terms 1 3^, 4(i for including the A isobar, the u>, rj, ao meson, and the 



cr [cf. Sec. 12.31 . The exchange potentials V are partial wave projected to the JLS -basis. 

The novelty in this work is the inclusion of the KY channels KA and K"L. This leads to a 
larger channel space and new transition potentials V to and within the KY channels. These new 
potentials V, related to the existing ones by SU(3) symmetry, are discussed in Sec. 12.31 They 
also contain the form factors which are used to regularize the scattering equation dj). 

2.2. s-channel processes 

In a model with explicit .s-channel states it is always possible to separate the amplitude into 
a pole and a non-pole part 

T = T p + T NP (2) 

where the pole part T p is defined as the set of diagrams that is 1 -particle reducible, i.e. there 
is at least one s-channel exchange. Usually, the non-pole, 1 -particle irreducible part T NP comes 
from t- and u— channel exchange processes collected into the non-pole potential V NP which is 
then unitarized using a dynamical equation of the type of Eq. dj]) — see also Eq. dj]) below. The 



separation of the type of Eq. (ffj) is widely used in the literature, see e.g. 142L 15411 . T NF is usually 
referred to as background, although the unitarization may lead to dynamically generated poles 
in r NP as discussed in detail in Ref. 15511 . There, the conclusion was drawn that the clearest 



separation into a background and a resonance part is given by the separation into a singularity- 
free part and the part fl_i/(z - zo) from the leading term in the Laurent expansion [cf. Eq. ©]. 

In the present study, we use the decomposition of Eq. (0, because the calculation of T F 
is numerically much faster than that of r NP . In a fit of only s-channel parameters, it is thus 
convenient to calculate T NP once and then fit the resonance parameters, which only requires 
the multiple re-evaluation of T v . Note that resonance w-channel exchanges contribute to all 
partial waves and are thus accounted to r NP . Nucleon, A, S, A(1232), and S*(1385) «-channel 



exchange diagrams are included with physically known coupling strengths [see Appendix A|, 
while H-channel diagrams from other baryonic resonances are neglected. Those would introduce 
additional parameters which are difficult to adjust for the diagrams do not introduce strong energy 
dependencies (for a discussion of u-channel contributions see Ref. J56ll ). 
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The pole contribution T p can be evaluated from the non-pole part T NP , i.e. from the set of 
diagrams that is 1 -particle irreducible. For this, we define the following quantities in a given 
partial wave, 

T NP (d,c) = V NP (d,c) + V NP (d,e)G(e)T NP (e,c) 
if (i, c) = c) + yf(i, d) G(d) T NP (d, c) 

r D (c,i) = j B {c,i) + T NP (c,d)G{d)y B (d,i) 
Z(iJ) = jf(i,c)G(c)r D (c,j) (3) 

where (Tp) are the dressed resonance creation (annihilation) vertices and £ is the self-energy. 
The indices i, j indicate the resonance in the case of multiple resonances, while c, d, e are indices 
in channel space. Integrals and sums over intermediate states are not explicitly shown in Eq. (O. 
For the two-resonance case, the pole part reads explicitly 15711 

r p = rD -i r m r = (r 1 ,r 2 ), rm = ( r 3, D = ( z " m I" 2u " 2l2 v \ (4) 

from which the one-resonance case follows immediately. The bare vertices y B for resonances 
with Spin / < 3/2 are derived from Lagrangians. The vertex functions for / > 5/2 are given in 
Eq. (IB . 3 1 > . Further details on the s-channel processes are given in Appendix B 



2.3. t— and u-channel exchange processes 

The t- and M-channel processes provide the non-resonant interaction in the meson exchange 
picture. The transition potentials without participation of KY have been derived in Refs. 13811391 
koll and explicit expressions can be found in these references. Here, we quote only the extension 
to the KY channels. The corresponding exchange processes are shown in Fig.Q] 

The vertices present in these diagrams are related to the already existing ones without strange 



particles using SU(3) symmetry (except for the <x meson, cf. Appendix A i. The coupling of 
SU(3) octets depends on two parameters which can be related to the axial coupling and an ad- 
ditional parameter. The values for these parameters have been taken from the literature and are 



not fitted in this study. This is explained in detail in | Appendix A There, one can also find the 



explicit amplitudes for the diagrams shown in Fig.Q] 

SU(3) symmetry is broken in the present study by the use of physical meson and baryon 
masses, as well as by different cut-offs in the form factors of the vertices. Exchange processes 
with strangeness S — -2 particles have been neglected because these baryons and, moreover, the 
corresponding 3-particle intermediate states, are heavy. A k exchange is in principle possible but 
not required by the data and thus has been neglected for simplicity. Furthermore, pN, 7tA <-> KY 
t- and M-channel transitions are neglected in the present work, as they appear only at loop order 
in the considered reactions nN — > nN and n + p — > K + l, + . 

2.4. Analytic continuation 

As argued in Ref. ll55ll . a clean separation of resonances and background is possible by the 
extraction of pole contributions from the analytic continuation. First results within different 



dynamical coupled channels models have been obtained in Refs. II41LI44L 15811 . 

The analytic continuation of the amplitude within the present framework has been derived 
in Ref. 14 111 in detail. Here, we summarize only the analytic structure. For the channels with 
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Figure 1: ;rA — » KY and AT — » KY transitions. For the other transitions used in this study, see Refs. 1391141 



stable two-body intermediate states, nN, r/N, KA, and KI,, there is one branch point at threshold 
^thres - m + ms which induces one new sheet. This is called the unphysical sheet. To search 
for poles on this sheet, it is convenient to rotate the right-hand, physical cut that extends from 
Zthres to co, into the negative Imz direction as shown in Fig. [2] Poles on this redefined sheet 
are close to the physical axis. Poles on other sheets are situated further away from the physical 
axis and thus typically have a much reduced effect on observables. However, there are certain 
situations in which such "shadow poles" can cause structures on the physical axis. An example 
is the Af*(1535). There, the interplay between the usual pole and the shadow pole causes the rjN 
cusp seen in the SI 1 partial wave, at least within the model of Ref. 114 111 . 

For the effective nnN propagators nA, crN, and pN, the analytic structure is more compli- 
cated: there is a branch point at z t h re s = 2m,, + which is induced by the cut of the self-energy 
of the unstable particle. Additionally, there are branch points in the complex plane at z t ' hres and 
tefhres)* w i*-h z thres - zo + M where M is the mass of the stable particle and zo is the pole position 
in the scattering problem of the unstable particle in the rest frame of the unstable particle 14111 . 
Those branch points can be regarded as pseudo-thresholds that have moved into the complex 
plane due to the unstable character of one of the particles. The argument to chose the direction of 
the cuts, associated with the branch points of the effective nnN channels, is the same as before: 
the cut is rotated into the negative Imz direction so that only those poles are found which are 
physically relevant. This is also indicated in Fig. [2] 

For the effective nnN channels, there are again situations, where poles on hidden sheets may 
have an effect on the physical axis: In Ref. 14111 a state in T NP in the Dn partial wave has been 
found, dynamically generated from the S -wave pN interaction (cf. also Refs. lf59L l60lo . This 
state, while its pole is well below the pN branch point, is visible as a washed-out structure at 
the nominal pN threshold, due to the fact that there is no direct connection from the pole to the 
physical axis, but only around the branch point at z' lhres in the complex plane. This resonance-like 
structure around z = 1700 MeV in the Do partial wave is, however, only visible in T NP ; once 
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Figure 2: Branch points of the coupled channels and the chosen directions of the associated cuts. Also, the isospin 
/ = 3/2 pole positions on the second sheet 7" <2 ' are shown [cf. Sec. 14. II . 



the strong N*(l520)Dn resonance is included in the same partial wave, resonance repulsion 115 511 
pushes the pole far into the complex plane filll . so that we cannot identify the dynamically 
generated structure with the three star A^*(1700)Di3 resonance 116 ill . 

Another example is the Roper resonance /y*(1440) whose poles are close to the nA pseudo- 
threshold in the complex plane. In case the coupling of the Roper to the ttA channel is large, this 
interplay of usual and hidden pole with the branch point may lead to the non-trivial structure of 
the Roper resonance visible on the physical axis 

The second sheet of the amplitude T denoted by J* 2) in the following is shown in Fig. |2 
with the cuts as defined above. In order to extract a pole residue on T (2) !_> -^ for a transition from 
channel i to /, we expand the amplitude J (2) in a Laurent series around the pole position, 



z-zo 



(5) 



In Appendix C the calculation of residues and branching ratios is discussed in detail, cf. Eqs. ( IC.6I ) 
and dC.8l ). respectively. 



2.5. Observables 

For elastic nN scattering, we compare with the dimensionless partial wave amplitudes r from 
the GWU/SAID analysis 113 ill . The scattering amplitude t for the transition ; — > / in channel 
space is connected to the amplitude T of Eq. (Q]) by 

kE to 

Tfi = -TT^PfPi Tfi, p= (6) 



where k (E, co) are the on-shell three-momentum (baryon, meson energies) of the initial or final 
meson-baryon system. 
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The observables in the reaction n + p — > K + l, + can be expressed via the t„n->ki.,i=3/2 ampli- 
tudes (abbreviated r in the following). According to, e.g. Ref. 16311 . for the scattering of a spin-0 
off a spin-^ particle, the differential cross section for the transition (k„ v) — > (kf, v') [v, V are the 
z-projection of the nucleon spin] can be expressed as 116311 

^ = \(s',v'\M\s,v}\ 2 (7) 
all 

whereas for the initial polarization 

Pi = <*vl<%v> (8) 

and for the final polarization 

p _ {M X M\MXv) (g) 
i {M Xv \M Xv ) ' 

Xv is the initial spin-| eigenvector and cr is the Pauli spin- vector. M can be written in terms of 
the non-spin-flip and spin-flip amplitudes g and h, 

M = g(k, 9)1 + h(k, 6)<r ■ h (10) 

where g and h are complex functions of the energy and scattering angle 0, and « = ^' x / . The 
polarization in the final state P j becomes 16311 

p = (\g\ 2 ~ \h\ 2 )Pj + (gh* + g*h + 2\h\ 2 ?j ■ n)h + i(gh* - g*h)Pj x h 
1 \g\ 2 + \h\°- + (gh*+g*h)n-P i 



while 



^ = [\g\ 2 + \h\ 2 + (g*h + gh*)h- P] -j-. (12) 



In the case of an unpolarized target, P, = 0, one obtains 



P (gh*+g*h) „_ 2Re(gr) , 

f \g\ 2 + W 2n \g\ 2 + \h\ 2n ( ^ 



and 



dcr t t kf 

-777 = i\g\ 2 + \h\ 2 )j-. (14) 
all kj 

The g and /i amplitudes can be expressed in terms of the partial wave amplitudes according 
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to 



1 

2yfkjki 



h = 



2 yjkjki 



Y{2J + 1) (d{ , (0) \t J{J -^- + T J{J+ ^-]sm - - d\_ , (6>) [r 
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The explicit expression for the differential cross section reads 
da 1 1 



da 



2k 2 2 

I 

1 1 

2£?2 



Y(2J+1)(t J(J -M +T J ^h-diA0) 

— ^ 7 3 



£(2/+l)(T 



)-d\ x (ff) 



and the total cross section is obtained by integrating over the solid angle Q, 

-4?Z (27+ 1)|t^' s '| 2 . 



1 JLS,L'S' 



(15) 



(16) 



(17) 



The spin-rotation parameter /3 is the rotation angle of the spin projection on the scattering 
plane. It is given by 164(1 



/} = arctan 



\\g\ -\h[ 



(18) 



3. Results 

3.1. Parameters and data base 

One bare s-channel state is included in each of the 1-3/2 partial waves S31, P31, D33, D35, 
F35, F37. Two are required by data in the P33 wave. These states were allowed to couple to all 
1-3/2 channels nN, K"E, nA and pN. Together with these four bare couplings, the bare mass 
has to be left free as a fit parameter. Thus, there are altogether 40 parameters for the pole part T v 
from Eq. ©. The values of these parameters can be found in Table |B.9l and the parameter errors 
are discussed in Sec. [5] 

The good description of the A(1232)f33 resonance shape requires also a fine-tuning of the 
cut-offs of the first s-channel state in P33, while for all other s-channel states, the cut-off was 
set to 2 GeV (cf. Appendix B 1. Additionally to the ^-channel parameters, the cut-offs of the 
diagrams of Fig. Q] were adapted (results may be found in Table IA.7l i. while those of the other t- 
and M-channel diagrams in the model ll40ll were not changed. 
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Figure 3: Contributions to the differential cross section for two typical energies: r NP (blue dashed lines), 7" p (magenta 
dash-dotted lines), and full solution (red solid lines). 



We fit to the n + p — > K + Y + differential cross section and polarization, given by the measure- 
ments of Candlin et al. II 1 Oil ,_available for z > 1822 MeV. For lower energies we have to resort 

. The latter data (see also Ref. |[7oTl ) are compatible 



67, 68 



to the data from Refs. J65L 

with the data of Ref. ft 1 Oil in the overlapping energy regions, but usually have larger errors. The 
polarization was re-measured later in the higher energy range of the considered data ll7lll . in 
consistency with the values of Ref. Iioll up to small deviations. See the captions of Figs. [4] to 
[9] for details. The spin-rotation parameter /? for the n + p — > K + I. + reaction has been measured 
in Ref. lfl5ll . Simultaneously, the polarization has been re-measured in Ref. Il5ll and consis- 
tency with results from Ref. ifioll was found. In summary, the considered data is consistent and 
represents the world data set from threshold to z — 2.35 GeV. 

For elastic kN scattering, the energy-dependent partial wave solution from Ref. 13111 up to F 
waves is used as input for the fit. Errors have been assigned to it by hand such that the nN data 
and the K + I. + data contribute similarly to the^ 2 . The uncertainties of the results presented in the 
following are discussed in Sec. [5] 



3.2. Differential cross section and polarization 

The differential cross sections for the reaction n + p — > K + I, + are shown in Figs. [3] S and 
The red solid lines show the result of this study. Overall, the data are well described over the 
entire energy range. For energies above 2 GeV, we do not claim validity of the present model, 
because the analysis of Refs. 11401 kill has been limited to that energy. Consequently, at the 
highest energies, the K + I. + data have not been fitted, but up to z ~ 2.25 GeV the description of 
the data is still good, as Fig.|5]shows. 

Note, however, that in the present Lagrangian-based framework, the amplitude allows for 
an extrapolation to higher energies; in analyses in which the potential is parameterized purely 
phenomenologically in terms of polynomials, there may be little control on the amplitude outside 
the fitted energy region. The overall agreement seen in Fig. [5] for the higher energies is good 
although a detailed inspection shows that there is room for some improvement. In particular, at 
energies > 2.2 GeV significant deviations are seen near cos 9 = ±0.5. 

To discuss the individual contributions to the differential cross section, we show in Fig. [3] 
for two typical energies, the non-pole part T NP (dashed line), T F (dash-dotted line) and the full 
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Figure 4: Differential cross section of n + p — » AT + E + from z = 1729 to z = 1985 MeV. (Red) solid lines: Present solution. 
Data: Ref. Q^l, except: z = 1729, 1757, 1789 MeV from Ref. |Hl,z = 1732, 1783, 1 8 1 3 MeV from Ref. (Hl.z = 1790 
MeV from Ref. (H, z = 1764 MeV from Ref. OH. 
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Figure 5: Differential cross section of n + p — > from z = 2019 to z = 2318 MeV. (Red) solid lines: Present solution. 

Data: Ref. [ljj. 
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Figure 6: Total cross section of n + p — > K + ~L + . (Red) thick solid line: Prediction from present solution. Also, the 
contributions from the individual partial waves are shown, as indicated. Data: filled triangles right from Ref. [10], other 
data: see references in the data compilation of Ref. [74}]. 



solution (solid line), see also Eq. In forward direction, the non-pole part T NP , i.e. the uni- 
tarized amplitude from t- and M-channel exchanges, produces a rise of the cross sections which 
becomes more pronounced as the energy increases and which is even stronger than the experi- 
mental forward peak at z — 2074 MeV. However, the resonance part T F produces a destructive 
interference with the T NP , which is crucial for reproducing the data. Note especially that T p is a 
lot more forward-backward symmetric than T NP at z = 2074 MeV. The forward peak shows the 
onset of the f-channel dominance which at energies > 3 GeV is most economically parameterized 
in terms of Regge exchanges 172l 17311 . 

In Fig. [6] the total cross section for the reaction n + p — > K + l, + is shown. The data have not 
been included in the fit, but the agreement is good. There is a slight underprediction of <x at 
z ~ 2 GeV by the present model, which comes from a slight underprediction of the forward peak 
in this energy range, also visible in Figs.|4]and|5] 

In Fig. [6] also the partial cross sections from the individual partial waves are shown. Except 
for the S31 and P33 partial waves, all other partial waves are very small; however, for the differ- 
ential cross section and polarization, their contributions is essential; indeed, while the removal of 
a resonance does not change much the total cross section, the differential observables can change 
drastically (see also Fig . [7j> . 

In Figs. [8] and |9] the polarization for the reaction n + p — > K + I, + is shown. Like for the 
differential cross section, the data show a rich and varying structure over the entire energy range, 
and the description by the present model is good. At energies z > 2.2 GeV, the data have not 
been included in the fit and are only plotted for comparison. We found the polarization to be 
especially sensitive to the resonance contributions, and the inclusion of these data is important to 
put constraints on the corresponding parameters. 

The influence of individual partial waves is illustrated in Fig.|7]for three typical energies. At 
lower energies, S —P wave interference is enough to describe the polarization, but not entirely the 
differential cross section (dash-dotted lines at z — 1813 MeV). For the latter, even a small F wave 
admixture is needed to explain the drop at cos 9 — -1 (cf. full solution). For z — 2019 MeV, one 
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Figure 7: Contribution from different partial waves to differential cross section and polarization for three typical energies. 
Dotted lines: S wave. Dash-dotted lines: S + P waves. Dashed lines: S + P + D waves. Solid lines: S + P + D + F + G37 
waves (full solution). 



also needs the D waves for an at least qualitative description of the polarization (dashed line), and 
the F wave is essential in the description of the details of the differential cross section. The same 
applies for the highest energy z = 2224 MeV, where all partial waves are needed to quantitatively 
describe the data. 

Fig. [TOl shows the spin-rotation parameter for the reaction n + p — > K + T. + [cf. Sec. 12.511 . In 
this study, ft is not included in the fit, but predicted (solid lines). The prediction from the isobar 
analysis of Ref. ITlll is also shown (blue dash-dotted lines). As ft is 2n cyclic, the data from 



Ref. Ill 111 (solid circles) have been plotted repeatedly (empty circles). The present model predicts 



P better than Ref. It 1 111 for z = 2107 MeV. Higher precision data would help further pin down 
the partial wave content because the results for /? already show that considering this observable 
is important to remove ambiguities in the partial wave content. 

3.3. Partial waves 

Figs. fTTI and [T2l show the / = 3/2 elastic nN — > nN partial wave amplitudes up to / = 7/2, 
except for G37 which is very small. 

The result of this study is indicated with the red solid lines. The data points represent the 



energy-independent partial wave solution from Ref. 13 1H . For comparison, also the previous 
solution from 2002 [40] within the framework of the Jiilich model is shown (green dash-dotted 
lines). Note that J > 3/2 resonances were not considered in Ref. lf40h . The contribution from 
r NP [cf. Eq. <J3J] is shown with the blue dashed lines. 



The description of the partial waves from the GWU/SAID analysis 13 111 is comparable to 



the results from the previous Jiilich analysis J40ll . While the present solution is better for, e.g., 
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Figure 8: Polarization of n + p — > A" + Z + from z = 1729 to z = 1985 MeV. (Red) solid lines: Present solution. Data: 
Ref. ful, except: z = 1729, 1757, 1789 MeV from Ref. (H, z = 1782, 1813 MeV from Ref. (H, z = 1790 MeV from 
Ref. (6J, z = 1764 MeV from Ref. (H, z = 1732 MeV from Ref. (H. 
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Figure 9: Polarization of n + p -> K + I, + from z = 2019 to z = 2318 MeV. Data: Ref. [T(J. (Red) solid lines: Present 
solution. 
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Figure 10: Spin rotation parameter ft of n + p — > K + l. + at z = 2021 and z = 2107 MeV. Note that ft is 2;r cyclic which 
leads to additional data points at shifted values shown by the empty circles. Data: Ref. fl31l . (Red) solid lines: Prediction 
from present solution. (Blue) dash-dotted lines: Prediction from Ref. flUl . 



the P33 partial wave, some deviations from the GWU/SAID analysis 13 111 at higher energies 
are visible in other partial waves. This may indicate the need for a more systematic fit of the 
parameters of r NP , or may be a sign of the tails of higher lying resonances. Note that similar 
problems for the elastic D35 partial wave have been in found in the GieBen i22ll analysis, and 
also in the EBAC analysis l7ol . 

Fig. Qj] shows the n + p — > K + I, + partial wave amplitudes obtained in this study (red solid 
lines). The contribution from the t- and M-channel processes, T NP from Eq. (f2]), are indicated 
with the (blue) dashed lines. The partial wave solution of Ref. fll ill is shown with the dash-dotted 
lines. Of course, the latter solution cannot be directly compared to the present one, because there 
is an overall undetermined phase. Still, even with such a global phase ambiguity, the figure 
shows that the partial waves are quite different. In particular, in Ref. Il 111 lower spin resonances 
are only included when providing a substantially improved^ 2 (cf. discussion in Sec. 14. lb . In the 
present analysis, we have used all those resonance states up to / = 7/2 needed to describe jiN 
scattering 13111 . Note, at least through coupled channel effects, they also couple to KH. 



4. Resonance analysis 

4.1. Pole positions 

In TableQ] the pole positions found in the present analysis are shown as Jiilich. The positions 
are visualized in Fig. [2] together with the chosen directions of the branch cuts (cf. Sec. 12.41 ). The 
first line of Table[T]indicates the data that have been taken into account in the different analyses, 
the second to fourth lines indicate the analyses (see below), their type, and whether the quoted 
values are pole positions or Breit-Wigner parameters. 

As for the well-established 4-star resonances, it is no surprise that the pole positions found 



in this study are in agreement with the values from the GWU/SAID analysis 113 111 (4th column), 
because the partial waves from that analysis serve as input for the present study. Indeed, the 
A(1232)F 33 , A(1700)D 33 , A(1905)F 35 , and A(1950)F 37 show clear signals in ttN nN (cf. 
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Figure 11: Elastic jtN — > nN partial waves S31, P31, P33, and D33. Data points: GWU/SAID partial wave analysis 
(single energy solution) from Ref. l3lll . (Red) solid lines: Present solution. (Blue) dashed lines: only T NP . (Green) 
dash-dotted lines: Jiilich model, solution 2002 from Ref. 1 40] . 
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Table 1: Pole positions zo of resonances with isospin / = 3/2. For each resonance, the upper row shows Re zo [MeV], 
the lower -2 Im zo [MeV]. The results of the present calculation are shown in the column Jiilich. The Table also specifies 
the type of analysis: Dynamical coupled-channels model [DCM], A^-matrix approach [KM], dispersion analysis [DA], or 
isobar analysis [IA]. The quoted values are either pole positions [P], Breit-Wigner value [BW], or pole positions obtained 
by speed plot techniques [SP]. For the other entries, see text. Uncertainty in the last digits in parentheses, (a) indicates 
that the corresponding resonance is dynamically generated in the present approach. 
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Figs. fTTI and fT2l and the present fit agrees well with these partial waves. The A(1620)5 , 3i is nar- 
rower than in the GWU/S AID analysis which comes from the weight of the KH data in the present 
fit. Also for the P3i(1910) resonance, the width is different from the one found in GWU/SAID; 
this may come from the small resonance signal in elastic nN scattering on top of a large back- 
ground (cf. Fig. [TTIi. However, note that the present formalism is not a /if-matrix approach; due to 
the dispersive parts present in this analysis, in principle one cannot expect similar pole positions, 
even if the amplitudes are very similar on the real, physical axis. 

The A( 1920)^33 and A(1930)Z)35 resonance show no or very small resonance signals in the 
GWU/SAID analysis of elastic nN scattering (cf. Table [2}. Their position is, thus, barely fixed 
from elastic nN scattering. It is then interesting to note that the constraints from the K + I. + data 
lead to resonance positions in vicinity to those quoted in the PDG, rated with 3 stars. Thus, we 
can accumulate further evidence for these states and their positions. It should be stressed again 
that the resonance positions are not preassigned in the present ansatz, but left completely free in 
the fit. 

Finally, we find poles in the scattering amplitude which are not induced by bare s-channel 
resonance states. Those poles are already present in T NP [cf. Eq. (fJJ] and arise from the uni- 
tarization of the t- and u-channel exchange diagrams. These dynamically generated poles are, 
in the present analysis, far in the complex plane: a A(1600)f33 and a A(1750)/ , 3i. Apart from 
these two poles listed in Table [1] we find a very wide dynamically generated pole in the S31 
partial wave at zo = 2170 - 645 ; MeV and one in the D35 partial wave at zo = 2734 - 445 i 
MeV which thus have widths F = -2Imzo of around 1 GeV. Both these states are too wide to 
be identified with resonances quoted in the PDG 16111 . This applies also to the P31 state, while 
for the A( 1600)^33 state quoted in Table Q] there may also be some evidence in the GWU/SAID 
analysis for a wide state. 

Note that not all those states included here by bare s-channel diagrams are necessarily gen- 
uine resonances; once the channel space is enlarged appropriately by inclusion of, e.g., res- 
onances like the A(1700)Z>33 may appear dynamically generated. This is discussed in Ref. lf75ll 
where the prediction of the I s and I c observables in the reaction yp — > n°j]p is shown to coincide 
well with experiment suggesting a dynamical nature for that resonance. 

The column of Ta ble [T] marked Giefien shows results of the /f-matrix based analysis from 
the Giefien group 1201 12211 . see also Introduction. While resonances with spin 5/2 have been 
included recently 12211 . the A( 1950)^37 resonance is not, which plays an important role in K + T. + 
production [11]. Also, the absence of some analytic properties restricts the model to real energies, 
and thus no pole positions can be quoted. The numbers shown in Table[T]are, thus, Breit-Wigner 
parameters. 

TableQ~]shows also the pole positions from the three standard partial wave analyses of elastic 
nN scattering, marked as GWU (George Washin gton University) Ir3lll . KH (Karlsruhe-Helsinki) 
lEQ, and CMB (Carnegie-Mellon-Berkeley) 112]. 

The followin g co lumn shows the pole positions from the extraction of the EB AC g roup as 
quoted in Refs. 144 16211 . based on the analysis of elastic nN scattering of Ref. 17611 . For a 
review on the theoretical foundations of the formalism, see Ref. 14211 . The framework has many 
similarities to the present one, although there are differences such as the treatment of the nucleon 
pole or the role of the Roper resonance, which appears dynamically generated in the present 
framework 13911 but is included as a genuine state in the EBAC model 06211 . 

The column DMT shows the recent pole extraction lf58ll from the Dubna-Mainz-Taipeh anal- 
ysis 14311 of elastic nN scattering. Like the EBAC and the present model, this approach is a 
dynamical meson exchange model, i.e., not a /if-matrix approach. 
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The column Cdl shows the Breit-Wigner parameters obtained from the isobar analysis of 



Ref. Illlll , In the isobar analysis, a purely phenomenological background and resonances are 
added in a way violating unitarity, and the fit is exclusively to the n + p — > K + Y, + data of Ref. fioll . 
The spin-rotation parameter ft has been evaluated in Ref. Jl5ll using this analysis resulting in 
poor agreement with the data (cf. discussion of Fig.|7]i. Relatively few resonances are quoted in 



Table[U because in the isobar analysis 11 ill only those lower-spin resonances are considered that 
lead to a significant improvement of the^ 2 . 

The column Mnly shows the results of the TT-matrix analysis of the nN — > nnN reaction of 



Ref. D24fl. 



The pol e positions from the multi-channel CMB type analysis of Vrana, Dytman, and Lee 



of Ref. 13411 are not shown in Table Q] This analysis finds two S31 resonances, two P31 and 
three P33 resonances, one D33 and one F37 resonances. It does not find the second D33 state of 
Table[T| but two D35 and two F35 resonances in contrast to the resonances quoted in TableQ] 

The last two columns show some of the predictions from the quark models of Loring, Metsch, 
Petry Q] and Capstick, Isgur Q]- In the work of Ref. 0], the non-strange constituent quark mass 
and two confinement parameters are fitted to the A-Regge trajectory. All mass values in Table [1] 
except the A(1232)P33 and the A(1950),F37 are then predictions. 

The question arises to which extent the poles found in this analysis can be related to the 
quark model states quoted in Table Q] J3, Si or others O, 77, 78]. As Table Q] shows, the quark 



models predict the baryonic resonance spectrum quite well. However, the hadronic dressing ef- 
fects are not explicit in these calculations, and they can be large. In fact, we have found large 
correlations between bare masses, coupling constants, the used channel space and the renormal- 



ization scheme 15511 . and the matching point between quark models on one side and dynamical 



coupled-channels approaches on the other side is still an open issue. 
4.2. Branching ratios 

In Table [2] the residues and resulting branching ratios into the nN channel are shown as 



obtained in the present study. The values are compared to the GWU/SAID results Bill [note 
that the values of r^/r tot from Ref. B~ltl have been obtained from a Breit-Wigner fit, while the 
values of the present study are directly obtained from the residues using Eq. dC.8H . 

As already noted in the discussion of the pole positions, most of the residue strengths and 
phases of the present analysis coincide quite well with those of the GWU/SAID analysis. As 
Tabled shows, this is the case for the 4-star resonances A(1232)f 3 3, A(1620)S 3 i, A(1700)D 33 , 
A(1905)/ r 35, and A(1950)F37. For resonances that couple weakly to the nN channel or which 
are very wide, such as the A(1910)P3i, the differences are larger. The dynamically generated 
A(1750)P3i quoted in Table[T]is too wide to be considered a resonance state. 

However, the dynamically generated A( 1600)^33 resonance is also seen in the GWU/SAID 



analysis Bill , with qualitatively similar properties (very wide, rather small to medium branching 
ratio into nN, similar residue phase 9). It should be stressed that the appearance of dynamically 
generated states in the present model is strongly restricted by the fact that the generating t- 
and M-channel processes connect all partial waves at the same time; there is little or no room 
to manipulate the strengths of these transitions in order to generate poles, without immediate 
consequences for all other partial waves. Thus, dynamically generated poles are rather stable 
objects — all of those found in this study are already present in the solution from 2002 of the 
Jiilich model l40ll . although they have not been searched for because the analytic continuation 
became available only in Ref. 114 111 . 
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Table 2: Left: nN — » nN residues of the present study (Ju) and from Ref. l3lll (GWU). For each resonance, the upper 
row shows \r\ [MeV], the lower 8 ["]. Right: nN branching ratios in %. (a) indicates that the corresponding resonance is 
dynamically generated in the present approach. 
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Table 3: Left: n + p — > K + T. + residues \r\ [MeV], 8 [°] of the present study. Right: Transition branching ratio [%] in the 
present study (Ju), from Ref. ITTll (Cdl), and from Refs. (20lE3l (Giefi). Uncertainty in the last digit in parentheses. 
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The attraction that leads to the dynamical generation of the A(1600)P 33 and the A(1750)P 3 i 
comes mainly from the nA channel: The coupling of the A(1600)P 33 into the (7rA)p 33 channel 
[cf. Eq. OJl] has a modulus of \g\ = 17 ■ 10~ 3 MeV~ 1/2 [A(1750)P 3I (ttA)p 31 : \g\ = 20] 
which is of the same size as the coupling of the A(1232)f 33 to nN (19 in these units). Indeed, 
the dominant decay channel of the A(1600)P 33 , quoted by the PDG JHH], is the one to nA (40- 
70%). Electromagnetic probes could shed further light on the nature of this resonance, as recently 
discussed using hadronic dressing and a constituent quark model 117911 . 

The A(1930)D 3 5 couples extremely weakly to the ttN channel in this analysis. Interestingly, 
the GWU/SAID analysis also finds a resonance in D35 with a very small signal in nN scatter- 
ing [cf. Table [2). Also the A(1920)P 33 resonance has a very small nN branching ratio in the 
present study — note that there is no pole found in the GWU/SAID analysis of elastic nN scat- 
tering. However, removing one or both of these states in the present analysis, the K + I. + data are 
described much worse, even if all other resonance parameters are refitted. 

In the second column of Table [3] the present results for the residues into the KH channel are 
shown. On the right-hand side, transition branching ratios of nN — > KI. from different analyses 
are displayed. Again, Ju marks the present results. 

1/2 1/2 

The values for (T nN FjL) /F tot are quite different from each other. The only common feature 
is the prominent role of the A(1920)f 33 resonance. Also, the F37 resonance has been found 



important in the present analysis and in Ref. Ill 111 



This wave is missing in the GieBen analysis which may distort the resonance content and 
branching ratios to KI, and may explain the differences observed in Table [3] As for Candlin's 



analysis II 111 , there are the above-mentioned conceptual problems of the isobar analysis, in par- 



ticular the oversimplified and unitarity violating construction of the partial wave amplitudes; 
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discrepancies to the present results are, thus, expected. 

There are also older analyses, not quoted in the Table [7(3, 8^1, which are based on low- 
statistic data previous to the one published in Ref. iflbTl . The branching ratios obtained in these 
analyses show large discrepancies among each other and also to those quoted in Table [3] 

All the discussed analyses differ in the data bases considered, the theoretical tools used, and 
third, by the quality of the fits. Moreover, the partial wave content is not unique even within 
the same framework and even if data from different reactions are combined into a global fit 
— for a clear demonstration of this, see Ref. 18 ill . Still, in the present approach, the explicit 
microscopical treatment of the non-pole part provides a realistic background which helps min- 
imize ambiguities from resonance contributions. Furthermore, the joint treatment of elastic iiN 
and K + Y, + data helps determine more precisely the resonance content of the K + l, + production 
amplitude. 



5. Uncertainties 

In this section we give some remarks on the reliability of the resonance parameters extracted 
based on the input data used — we will make no attempt to estimate the theoretical uncertainty 
of the approach as such. We are not (yet) able to quantify the uncertainty introduced into the 
analysis by the particular formalism used. In principle, once a set of model analyses exists fitted 
to the same data with the same channels included but based on different formalisms, a comparison 
of the resonance parameters extracted should provide this information. 

As mentioned in Sec. 13.11 the error bars in the x 2 minimization have been taken from exper- 
iment for the reaction n + p — » K + H + , but assigned by hand for the partial waves of elastic jiN 
scattering, since no uncertainties are provided for the energy dependent partial wave amplitudes 



provided by the GWU/SAID analysis 13 lH and the uncertainties provided for the corresponding 
energy independent analysis do not have direct statistical meaning. The uncertainties for the nN 
partial waves are chosen such that the contributions from both reactions to the total x 1 are ap- 
proximately equal. This makes a rigorous error analysis of the present results impossible, which 
would require a fit directly to the elastic nN data. Nevertheless, assuming these assigned errors 
are realistic, we outline in this section how to obtain in principle the uncertainties on the parame- 
ters and derived quantities, like pole positions and residues. The error (0.01) and energy spacing 
(40 MeV) used to include the nN partial waves in the^ 2 minimization are shown in Fig.ll4lfor 
the example of the F35 partial wave. 

As mentioned in Sec. 13. II 40 parameters tied to the resonances have been varied to minimize 
the x 2 - There are also other parameters tied to the non-resonant part, given by the form factors 
shown in Table IA.7I The latter have been roughly adjusted by hand before carrying out the 
numerical fit of the resonance parameters, and we do not consider them as free parameters for 
the error analysis carried out in the following. In the space of 40 parameters the error of parameter 
Pi is determined by the range of p, for which the best^- 2 . rises by less than Ay 2 = 1, optimizing 
at the same time all other 39 parameters. In the limit Ay 1 — > 0, this non-linear parameter error 
approaches the usual parabolic error that can be obtained, e.g., from the Hesse matrix. 

Here, we consider only the example of the F35 partial wave. We restrict the determination of 
the error to the 5-parameter subspace tied to the F35 resonance, e.g., its bare mass [cf. Eq. (|4]i] 
and four couplings to the channels nN, pN, nA, and KI. [cf. Eq. (f3]l, | Appendix B| . Furthermore, 
we determine the parameter errors within this subspace. This means that in the optimization 
involved in the determination of the parameter error (see above), only the 4 parameters of the 
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Figure 14: Left: The F35 partial wave in nN — » nN. Data points: error given to the energy dependent SAID solution f3lll , 
as used for the numerical fit. Thick (red) line: Minimal x~ solution. Thin (black) lines: Representative solutions in the 
X 1 + 1 criterion (determination of the non-linear error). Center and right: The A(1905)F35 pole positions and branching 
ratios from those solutions. The ellipses are introduced to guide the eye. 



Table 4: Error estimates of bare mass m\, and bare couplings / for the A(1905)F35 resonance. For other bare parameters 
see TableDT9l 



m b [MeV] 



2258!$ 



7TN 

0500+ 00011 



pN 
-1 62 +L29 

1-u -1.61 



nA 

i i c +0.030 
i- -0.022 



120 +00065 



subspace are varied while leaving the other 35 at the optimum. We have checked that this restric- 
tion has surprisingly little influence on the parameter error because partial waves are explicitly 
included in the fit instead of nN observables; for example, varying a resonance parameter of the 
P33 partial wave influences a parameter error of the F35 partial wave only indirectly through the 
inclusion of the KI. data in the total x 1 - 

The non-linear parameter errors are shown in Table [4] The errors are small for the bare 
couplings to the nN and KH channels, for both of which the data constrain the values. The errors 
are larger for the bare coupling to the pN state that is less constrained by data. Indeed, there is a 
strong correlation between the bare pN coupling and the bare mass, which therefore also has quite 
a large parameter error. Although no data are included for the nA channel, the corresponding bare 
coupling has small errors. This is because the nA channel provides most of the nnN phase space 
that is responsible for the inelastic resonance width, which is well constrained by the elastic nN 
amplitude as shown in Fig.[14]to the left. 

From the non-linear parameter errors, one can determine the uncertainties of derived quan- 
tities such as pole positions, residues, branching ratios, or the amplitude itself. To scan the 
parameter space within the errors, we have taken four sample points of a given parameter within 
its error (always optimizing all other parameters). As there are five parameters in the considered 
subspace, 20 solutions are obtained from which the F35 amplitude and the A( 1 905)7^5 pole po- 
sition and residues have been calculated [see Fig. [141 and Table |5). The error on these quantities 
is then given by the maximal range reached by these solutions. In Fig. [141 we show also that 
there are correlations between real and imaginary part of the pole position and also between the 
branching ratios. As mentioned before, the transition branching ratio is indeed better determined 
than the individual branching ratio into K"L; we obtain (T 1 ^ 2 r^)/F tot = %. 

Note that the F35 amplitudes allowed by the discussed^ 2 + 1 criterion, shown in Fig.[14]to the 
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Table 5: Error estimates of pole position and residues for the A(1905)f35 resonance. 







nN -> nN 


nN -> fCZ 


Re zo [MeV] 1764+^ 


\r\ [MeV] 


U + 1.7 
11 -1.4 


1 4+0.24 
1 -0.21 


Imz [MeV] -l09 + _\ 3 2 




-45+3-8 

-11 


-313+ 4 - 2 



left, lead to a much larger rise A¥^F35 ^ 1 m me A^\rF35 °^ tne n ^ c ' ata a l° ne ' as an inspection 
by eye shows. Still, these solutions fulfill the x 1 + 1 criterion due to the contribution from the 
KH data to the total x 1 - If one determines, e.g., the uncertainty of the pole position or branching 
ratio from nN data alone, one would obtain, of course, much smaller errors on these quantities. 
Thus, the uncertainties on pole position and branching ratio shown in Fig.[14]and Table [5] should 
be understood as upper limits. 

In summary, we have outlined how to determine the statistical errors of the present results, for 
the example of the F35 resonance. A rigorous statistical analysis, as outlined above and carried 
out in Ref. [82] for K p scattering, requires a direct fit to nN observables and the full inclusion 
of KA, KI,(I =1/2) and r/N data and will be carried out in the future, but the present discussion 
serves to illustrate the error one expects from such an analysis. 



6. Summary 

A first combined analysis of the reactions nN — > nN and n + p — > K + l, + within the unitary 
dynamical coupled-channels framework has been presented. For the n + p — > K + l, + reaction, the 
world data set from threshold to z = 2.3 GeV has been considered. 

Dynamical coupled-channels models are particularly suited for combined data analyses: the 
SU(3) flavor symmetry for the exchange processes allows to relate different final states. The t- 
and M-channel diagrams connect also different partial waves and the respective backgrounds. 

As a result, for both nN and KL, a realistic and structured background can be provided, de- 
pending only on a few free constants and form factors whose values are all in a natural range. 
Consequently, only a minimal set of bare ^-channel resonances is needed to obtain a good fit to 
the combined data sets. This may also be tied to the fact that in this field-theoretical, Lagrangian 
based approach, the dispersive parts from intermediate states are fully included and thus, analyt- 
icity is ensured. 

Apart from the well-established 4-star resonances, a wide A( 1600)^33 state has been found, 
dynamically generated from the unitarization of the t- and M-channel exchanges. Furthermore, 
there is a clear need for the three-star A(1920)P33 resonance. This state is found to couple only 
weakly to nN but stronger to KH. Thus, in the present combined analysis of elastic nN scattering 
and K + T, + production, evidence for a "missing resonance state" |Q] could be accumulated which 
indeed has no clear signal in elastic nN scattering alone. 
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Figure A. 15: t- and w-channel exchange processes. 
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Appendix A. Exchange potentials with KY states 



In |Appendix A.l we list the explicit expressions for the t- and u-channel exchange diagrams 



that involve the KA and KL channels. For the other exchange processes contained in the model, 
see Ref. l40ll . The new coupling constants for the exchange processes with KY participation are 
related to the cases without KY through SU(3) symmetry. The corresponding expressions can be 



found in Appendix A. 2 



The kinematical quantities are specified in Fig. IA.15l The index 1 and 3 (2 and 4) denote the 
incoming and outgoing baryon (meson). The on-shell energies are 



for the baryon and the meson, respectively. In the TOPT framework used in this study, the zeroth 
component of the initial and final momenta are set to their on-mass-shell values: p9 = £, or 

p° = CJj. 

q is the three-momentum of the intermediate particle, q with q° = E q (baryon exchange) or 
q° = u q (meson exchange) means the 4-momentum in the first time ordering whereas q indicates 
the second time ordering with q a = -E q (baryon exchange) or q° = -cj q (meson exchange). 
Furthermore, in the potentials quoted in | Appendix A.l| F™ is the Rarita-Schwinger propagator 
of spin 3/2 particles given in Ref. ll39Tl . and p = y 11 p^. 
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Table A. 6: Isospin factors for exchange diagrams with participation of the KY channels. See Refs. [39. 40] for the 
corresponding values for the other diagrams. 



Transition 


Exchange 


IF(i) 




Transition 


Exchange 






ttN -> AK 


K* ex. 


V3 





AK -> AK 


cr ex. 


1 







X ex. 


VI 







ai ex. 


1 







X* ex. 


V3 







(p ex. 


1 





nN -> T.K 


K* ex. 


1 


2 


AK -> ZK 


X"* ex. 


-V3 







X ex. 


2 


1 


2AT-> Xtf 


cr ex. 


1 


1 




A ex. 


-1 


1 




w ex. 


1 


1 




X* ex. 


2 


1 




ex. 


1 


1 












p ex. 


2 


-1 



If the intermediate particle is a X* baryon, we use 118311 q^_ t — e\ — €4 with 



999 999 
Z + m, - mi z-mi + mi 

ei = T -• 64 = T ~ (A - 2) 

2z 2z 

Each exchange diagram includes a kinematical normalization factor 

K=-r^-r — (A.3) 

(2ny 2 yZIi^M 

The isospin factors IF for the exchange processes without the participation of the KY channels 
can be found in Refs. j3^ . l40ll . For the extension to the KY channels, carried out in the present 
work, the isospin factors can be found in Table IA.6I Also, every exchange process quoted 
in Appendix A.l is multiplied with form factors corresponding to the two vertices, quoted in 



Appendix A. 2 



Appendix A.l. Amplitudes for the exchange diagrams with KY 

The following expressions give the amplitudes shown in Fig.[T]in the helicity base, 

(A'\V(k,k',z)\A) (A.4) 

with the notation for the momenta as in Fig. lA.15l and the dependence of the amplitude is on the 
cm. (off-shell) three-momentum of the incoming (outgoing) meson-baryon system, k(k' ), and 
the total scattering energy z. To solve the scattering equation in the JLS -basis, these expressions 
still have to be partial wave projected |38, 39[ 4oh . 
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ttN -> KA 



• /T* f-exchange 



KgKnK- u{_p 3 ,A 3 )(\\g KNK ,"f ~ ; ™ K O*™ q v \ 



2lflN Z — Oi q — £3 - W2 

+ [ghNK'J^-i- 0* v q v ] )u{puh) IF 

2niN z - (x>q - E\ - ti)\> 2u) q 



I. M-exchange 



gNiKgniA ^ , ,.,5 q + mz q + mz 



k 2 u(p 3 ,A 3 )y — - + - - —\y p 4 u(p l ,A l )IF 

mj, 2E q \z -E q -co 2 -o)4 z- E q - E\ - E 3 / 

• S* M-exchange 

gm.>KgnL>A . . Plf^q-L*) ( 1 1 

< 2 u (P3'h) ^ p + p p ^r\P4yU(puM)IF 

mi 2E q \z- E q - cd 2 - cd a z-E q -E x - E 3 ) 

nN -> KZ 

• K* r-exchange 



KgKnK' U(p 3 , A 3 )([gz N K>y ~ i-^-O^qv] 



2m N z - (O q - £3 - u>2 
+ [gLMr-7^ - ^ ?v] )u(p u A x ) — IF 

2l1lN Z — Oiq — E\ -0)4' 20) q 



2 M-exchange 



* 5 «(P3./l3)r ^ p + s p 7 p4u{p\,A\)IF 

m l n 2E q \z - E q - 0)2 - w 4 z-E q -E\-E 3 



A M-exchange 



gNAKgnAY. -,->,. 5 P2 q + m A q + OT A 
C 2 «(P3, ^3)7 TTT" p + p p 7T 

m£ 2£ 9 \z - £ 9 - a> 2 - m Z- E q -E\- E 3 



.5 , 



• Z* M-exchange 

K 2 M (P3>40 p + T, T, —)P4 r U{Pl,M)IF 

mi 2E„ \z- E„- ah -04 z-E q -E x -E 3 
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KA -> KA 

• cr f-exchange 



ga-AAga-KK 

<- — : ; t. T. + T. T. ]u(p 3 ,A 3 )u(pi,A 1 )IF 

2m n 2co q \z-E q -E3-(02 Z- E q - E\ - 0)4 

• <d f-exchange 



KgcoKK U(P3, ^([gAA^ ~ jy^^tfv] 



2m N z - u) q - £3 - o> 2 

+ [gAAwT-ii: q v ] IF 

2mti z - Oiq — E\ - L04' 2<x> q 



• f-exchange 



Kg<f>KK u(Pl, ^([gAA^ ~ i^^-^qy] l — 

v 2m N Z - Oiq - E3 - 0)2 

„ JAA4, 1 \ rf (P2 + P4)» 

/nixj 7 — — h.\ — ill* ' 



IF 



2mjv z — tOq - E\ —0)4' 2u) q 

KA ^ KZ 
• p f-exchange 



K gKpK U(p3, /^([gApx/ - i^^Cr^ V ?v] ~^ 

v 2m N z-OJq- £3 - w 2 

+ [^ApS-r-'T; <r q v \ )u(p u Ai) — IF 

2niN z- U) q - h\ - u)\' 2oj q 



KL -> KL 



cr f-exchange 



gKo-KgZaZ (- 2 P2,f4) / 1 1 \ 

k — o ~ ^ + i i m(p3,^3)m(Pi,^i)/£ 

2f7Z,r 2w ? \Z-Eq-E3- 0)2 Z~ E q - E\- Oii,) 

• w f-exchange 



2m N z-(Oq- E 3 - C02 

+ [giZcoT-i-z q v ] )u(pi,A 1 ) — IF 

2mjy z — oi q — E\ -0)4' 2u>q 
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• 4> f-exchange 



Kg^KK u(p3, ^)([gi£4,y - i^-a^qy] ^ 

v 2m N z-(O q - E3 - (jJo 

+ [gizfT-i-z V* q v ] )u(pi,M) z IF 

2niN z - a>q - 1 1 -0)4' 2a> q 



p f-exchange 



fXLp - 



1 



+ [g^ P y - h — ^ v qv] 

2m m 



a>4 



- (D q - £3 - U 2 

\,r* 1 ^ 0*2+ Mi 

jwCPi.'ti) — = IF 

' 2u n 



Appendix A.2. Coupling constants, form factors (t-, u-exchanges) 

First, the coupling constants needed for the exchange diagrams (Fig.Q] expressions in Appendix A. 1 
are quoted, followed by the form factors and their values. 

The coupling constants are related to their counterparts without strange particles by SU(3) 
flavor symmetry l84ll as outlined in Sec. 12.31 The coupling 8b ® 85 ® 8m of baryon octets 8b and 
pseudoscalar meson octet 8m depends on two parameters g\, g2 corresponding to the symmetric 
and antisymmetric octet. They can be related to the coupling constant g and the mixing parameter 
a in the notation of Ref. [84], which is also used here, 



V30 V6 



40 



fi+ 24*2, 



V6 g2 

24 7 



gi and g2 can be also expressed in terms of the standard D and F couplings [61]: 



D 



V30 



V6 
24 ; 



(A.5) 



(A.6) 



with a = F/(D + F) (definition of a of Ref. 18411 ). The couplings of the physical u> and (f> are 
obtained assuming ideal mixing of the SU(3) states lo\, oj%, i.e., the <p meson does not couple to 
the nucleon (see, e.g., Ref. Il6ll0 . 

Within SU(3), the mixing angle a is a free parameter but can be determined from SU(6). We 
use Il85ll (The index P denotes a pseudoscalar meson and V means a vector meson) 

cibbp = 0.4 , a BBV = 1.15 , a PPV = 1 . (A.7) 

For the coupling of vector mesons to octet baryons we use the following Lagrangian: 

Lint = -gNNpViy ~ ^^ V d v ]?p M (A.8) 

2m^ 

which consists of a vector part with y and a tensor part with o 4iv . The three fields *P, ¥ and p 
are connected through a vector coupling g NNp and a tensor coupling f BBV = gNNpKp US 87]. The 
SU(3) couplings quoted below are therefore divided into a vector part and a tensor part. We use 
k p = 6.1 li!] and f ma = 0. 
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The coupling of the A and the S to the cr meson was not fixed by SU(3) symmetry since 
the cr is not a stable particle but regarded as an effective meson-meson state (38, 39[ 4(i, i.e. as 
a correlated nn exchan ge proc ess. In the calculation for the background diagrams we take the 
following values from 149118611 . obtained from the hyperon nucleon interaction: 



gAAtr = 8.175 , gKKcr = 1.336. 



The coupling g^za was set to a different value compared to 149118611 . We choose gzz 
The SU(3) symmetric coupling constants are given by the following expressions: 

• Couplings for octet baryon, octet baryon and pseudoscalar meson: 

gNNn 
gUNK 



(A.9) 
29.657. 



= gBBp(i -2a BBP ) 



ghNK 

gM.1T 



\ V3g BB p(l +2a 



BBPI 



V3 



gBBp(l ~ CIBBP) 



• Vector coupling for octet baryon, octet baryon and vector meson: 

gNNp = gBBV , 

gNNoi = gBBv(4aBBV ~ 1) > 

gANK' = -gBBV^i^ + 2aBBv) > 

gl.NK' — gBBvi^ - l-OCBBv) > 

gAAoj = gBBV^(5aBBV ~ 2) , 

gU.0) = gBBv2&BBV i 

gAA<p = -gBBV-Y(2-ClBBV + 1) > 

gY£4 ~ -gBBV^2{2a BBV - 1) , 

gZZp — gBBV^aBBV , 

gAZp = - &BBV) ■ 

• Tensor coupling for octet baryon, octet baryon and vector meson: 



Inn p - 


gNNpKp , 


fANK' = 


f 1 f V3 
-jNNui - JNNp — 


fzNK> = 


—fNNai 5 + fNNp 2 ' 


/aAoi = 


fNNdf, - fNNp J > 


JyLlo = 


fNNoJ J + fNNp 5 i 


/aa* = 


~fsN<o Jy= - /iVJVp 


faty = 




fel-p = 




fAUp = 


"JWiViJ ^ + JNNp — 
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(A. 10) 



(A. 11) 



(A. 12) 



Coupling for pseudoscalar meson, pseudoscalar meson and vector meson: 



gimp = 2gppv , 

gKKp = gPPV , 

gKnK' = -gPPV i 

gKKai = gPPV i 

gKKtp = ^flgPPV ■ (A. 13) 

• Coupling for decuplet baryon, octet baryon and pseudoscalar meson: 

gANx = gDBP , 
gZ'NK = ~gDBP-^ , 

gZ'ln = gDBP-^ , 

gZ-An = gDBP^ ■ (A. 14) 

For the new diagrams of Fig. Q] we use the same expressions for the form factors as in 
Ref. H, 

A 7 1 " 

A z - mr r 
A 2 + q 



F (q) = \ A , , J ) (A. 15) 



where m x is the mass and q the momentum of the exchanged particle, n — 1 , 2 denotes a 
monopole or dipole form factor. The dipole type applies to vertices with p, A, K* or 2* as 
exchanged particle, otherwise a monopole form factor is used. For the form factors of all other 
exchange diagrams of the model, see Ref. H|. The numerical values for the new form factors 
are given in Table I A. 71 The numerical values for the other diagrams (without KY) have not been 
changed, see Ref. 14011 for the values. 



Appendix B. Bare resonance vertices 



The bare resonance vertices for J < 3/2 are given by the effective Lagrangians listed in 
Table IB. 8l There, p^ v = d^py - dvPp for the field of the p meson. The vertices derived from these 
Lagrangians are partial wave projected to the JLS -basis 1 39[ 4oTl . 

The Lagrangians for the KA couplings have the same structure as for the r]N couplings, the 
ones for the K"E couplings have the same structure like for the nN case except for the replace- 
ments T7? — > tTj; , sft — > or 5*7? — > Thus, they are not quoted explicitly in 
TableES 

Additionally, form factors are supplied for each vertex, given by 



F(k) - 



A 4 + m\ 



A 4 + (E(k) + u)(k)) A j 



(B.l) 



where m# is the nominal mass of the resonance (~ Rezo) and E(k), a>(k) denote the on-shell 
energies of the incoming or outgoing baryon and meson with cm. (off-shell) momentum k. We 
have n - I'm the case of J < 3/2 for all channels except for A7r (« = 2). For J > 5/2 we have 
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Table A. 7: Form factors A for the exchange diagrams with KY. The columns Ex specify the exchanged particle. For the 
numerical values of the form factors of the other diagrams, see Ref. (ioll . 



Vertex 


Ex 


A IMeVl 


Vertex 


Ex 


A rMeVl 


ttK*K 


K* 


1700 


KcrK 


O" 


1400 


NK*A 


K* 


1200 


KcoK 




1600 


NK"Z 


K' 


1800 


K(pK 


<f> 


1500 


ttZA 


I 


1800 


Ao-A 


(T 


1000 


ttI*A 


r 


2000 


AcoA 


CO 


2000 


nTL 


s 


1800 


A<f>A 


<P 


1500 




s* 


2000 


ApE 


P 


1160 


NXK 


z 


1800 


ZcrZ 


(T 


1000 


NI>*K 


z* 


2000 


ZwZ 


OJ 


2000 


ttAE 


A 


1800 


202 


</> 


1600 


NAK 


A 


1800 


ZpE 


P 


1350 



n — 2 for all channels except for An (n = 3). For the cut-off parameter we choose A = 2 GeV, 
except for the A( 1232)^33 resonance, where we fine-tune A = 1.8, 1.7 GeV for the vertices to 
the 7rAf, ttA states, respectively. 



In Appendix B.l we list the partial wave projected, bare resonance annihilation vertices v. 
The bare coupling constants obtained in the fit can be found in Table IB. 9l Every vertex function 
is multiplied with the corresponding form factor of Eq. ( IB. lb and the isospin factor Ir listed in 
TableETO] 



y B = F(k)I R vJ^p^- (B.2) 

with jb from Eq. ([3]). The resonance creation vertices y-g are given by = (y^) ' . In Eq. dB.2b 
and Appendix B.l| E and to denote the baryon and meson on-shell cm. energies, respectively, 



niB is the baryon mass of the channel, and Ic = \k\ is the baryon-meson cm. momentum. 

As mentioned in Sec. 12.21 the resonance vertices with total spin J > 5/2 have not been 
derived from Lagrangians. Instead, they have been constructed obeying the correct dependence 
on the orbital angular momentum L (centrifugal barrier). From parity considerations, one can 
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Table B.8: Effective Lagrangians for the resonance vertices. 



Vertex 


■Lint 


Vertex 


■Lint 


N*(S i 


i)Nn 


^N-yftd^v + h.c. 


A*(S 3l )Nn 


-LAyf^djy + h.c. 

win " 


AT (Si 


i)Nri 


i^N'T^d^ + h.c. 


&*(S 31 )Np 


/A*SV rV/^ + h.c. 


N*(S 


ONp 


f^N-y'y^P^ + h.c. 


A*(S 3 i)An 


■J-A*y 5 fdjAV + h.c. 


N*(S 


i)Aw 


^•y^c^A" + h.c. 


A*(P 31 )Nn 


-^AVy^^T + h.c. 


N*(Pi 


i)Nn 


-£^-rV^^7?^+ h.c. 


A*(P 31 )N P 


-fA*^'y^p ll x ¥+ h.c. 


N*(Pi 


i)Nr, 


-£ < *VrVc^ v i' + h.c. 


A*(P 3l )An 


^A*fdjA^ + h.c. 


N*(Pi 


ONp 




A*(P 33 )Nn 


£-A*&d"7t*¥ + h.c, 


N*(P 


i)Att 


^^d^A" + h.c. 


A*(P 33 )Np 


-^A^Vv^'p^^ + h.c. 


N*(Pi 


3 )Nn 


£%,fd M 7tV + h.C. 


A*(P 33 )An 


^AtffTdjAV + h.c. 


N*(Pi 


3 )Nt, 


£^.C^ + h.C. 


A*(D 33 )Nn 


4A;,r 5 r,i t Wif , p+ h.c. 


N*(Pi 


3)Np 


-^^•rVW? + h.c. 


A*(D 33 )Att 


i^f/^A' + h.c. 


N*(P 


3 )An 


£^, r V^3 v 7?A p + h.c. 


A*(D 33 )Np 


^A; rv 5 / V iV ^ + h.c. 


N*(D 


3 )Nn 


±Vy 5 r?d v d^ N , + h.c. 






N*(D 


3 )Nr] 


^y 5 y v d v d,rj%, + h.c. 






N*(D 


3 )Np 


£*£.r v WP + h.c. 






N*(D 


3 )An 


i^%Jrd^A v + h.c. 







Table B.9: Bare resonance parameters: masses and coupling constants /. 





m b [MeV] 


nN 


pN 


7TA 


IX 


A(1232)P 33 


1535 


1.44 


5.88 


-0.551 


0.0316 


A(1620)5 3 i 


3669 


0.769 


1.107 


-6.05 


2.25 


A(1700)A,3 


3442 


0.100 


-6.47 


-0.845 


0.170 


A(1905)F 35 


2258 


0.0500 


-1.62 


-1.15 


0.120 


A(1910)P 3 i 


3114 


0.367 


4.36 


-0.355 


0.231 


A(1920)P 33 


2508 


-0.123 


-2.96 


-0.530 


-1.86 


A(1930)D 35 


2332 


0.177 


-4.19 


-0.178 


4.12 


A( 1950)F 37 


2597 


0.580 


12.3 


1.87 


0.663 
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Table B.10: Isospin factors Ir for resonances vertices. 



Nn 


Nrj 


Np 


An 




T.K 


I=\ V3 


1 


V3 


-V2 


1 


-V3 


'=1 1 





1 







-1 



easily relate the bare vertices from the resonances with J < 3/2 to those of higher spin, 

U, \MB _ ^ / \MB /• \MB _ ^ / \MB 

2 Mb 2 2 m B 2 

U, \MB _ ^ / , .MB Ul nMB _ ^ , nMB , r ^ 

2 fl!j 2 2 OT^ 2 

with the {jE^f? from Eq. dB.2l ). Eq. ( IB. 31 ) provides the correct dependence on L for all channels 
MB = nN, pN, nA, and FL. 

Appendix B.l. Partial wave projected resonance vertex functions 

In the following, the classification of vertices corresponds to the quantum numbers of the 
nN channel. E.g., 511(531) refers to the I - J - 1/2 resonances that couple to nN in 5-wave, 
as for example the iV*(1535)5 n or A(1620)5 , 3i. The other channels can, of course, couple with 
different orbital momentum L and S (note there are three pN and two nA channels with different 
combinations of L and S). 

In all cases, the vertex functions v for the resonance couplings to KY and r]N are the same as 
for the nN case except for the different masses and isospin coefficients from Table IB. 101 Thus, 
the r]N, KY vertices are not quoted explicitly. 

Sn (S31) 

• Nn 

f 1 

v = - i ■=- (u n +E N - m N ) 

mn V8 7T 

• Np (L=0, S=l/2) 

v = - / — — (cj p + E N - m N + 2 m p J 

ylAn m p 

• Np (L=2, S=3/2) 

u = - f 

y[\2n m P 

An 

. f 1 E h + u n 

v - -1 — (E A - m A ) 

m n ^]Y2n m A 



11/ X 

v = - / — — \a> p + E N -m N - m p ) 
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Pu (P31) 



Nn 



Np (L=l, S=l/2) 



iVp (L=l, S=3/2) 



v = i — U 



m„ a/8 7t I E N +m Nl 



yf2ATt m p \ E N + m N E N + m N l 



2 m n 



1 k 



1 + + 



VT2 7r m p \ e n + m N E N + m N , 



An 



P 13 (P 33) 

• Mr 



Mo(L=l,S=l/2) 



iVp (L=l, S=3/2) 



/ 1 k 
v - 1 — (E A + u>„) 



v = — 1 £ 



^ V24tt 



v = 



/ 1 



(E N + m N - a>p - m p ) 

+ m N v hp/ 



v - 



/ 1 



Np (L=3, S=3/2) 



An (L=l) 



"V V3607T E N +m N 



f 1 



^5 £V + 5 niN + 4 Dp + m p ) 



7= t; — ; (w p - nip) 

/ 1 



An (L=3) 



m * V360^ '"a 



/ 1 k 



(£ A + 4m A ) 1 + 



E A +m A i 



v = -i — (E A - m A ) 1 + 

m n V407r m A \ E A +m A) 
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Dl3 (P33) 

• Ntt 

• Np (L=0, S=3/2) 

• Np (L=2, S=l/2) 

• (L=2, S=3/2) 

• Att(L=0) 



/ 1 

v = — 7 . — (E N - m N ) (to„ + E n + m N ) 
K V247T 



. / 1 

v = - 1 



V72tt 



. / 1 

v = - 1 



m p V72tt 



v = — i — ; (<jj,, - m + 2En — 2111^] 

mp -fan ' 



(oj p - m p - E N + m N J 



f 1 1 

v = — (a^ + E A - wa) (-Ea + 2 m A ) 



• An(L=2) 

f 1 1 

v = - — (oj„ + E A - m A ) (E A - m A ) 

m n Sl2ii m A 

Appendix C. Residues and branching ratios 

The residue a_i and constant term ao from the Laurent expansion of Eq. (|5]l can be obtained 
by a closed contour integration along a path F(z) around the pole position zo, 

1 r TV\ Z )dz fru 



Alternatively, the residue and subsequent terms in the Laurent expansion can be obtained by an 
iterative procedure according to 

= - 

dz" r®( Z ) a-i 

d 1 , 1 2fl 



3 3 1 6(«q — 



(C2) 



which is numerically stable (the inverse T (2) matrix has a simple zero at z — Zo) and fast (no 
integration required). 
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The residue and constant term, a_i and ao, can be expressed in terms of dressed quanti- 
ties Q, 



■it) 



ao - 



r r (t) 



(|(r D r^ ) + ^^z) + r^W. (c.3) 



where To (r^', £) is the dressed annihilation vertex (creation vertex, self-energy) as defined in 
Eq. (0, evaluated on the second sheet at zo- Eq. ( 1C.3I ) shows that there is a contribution to the 
constant term aq from r NP , as expected, but also from the pole term. This is one of the reasons 
why an identification of T NF as background is problematic [55]. 

For the two-resonance case shown in Eq. ©, the residues can be expressed in terms of 
dressed quantities like in the one-resonance case of Eq. iC.3i . 



a-U = 



det£> (2) 
d/dz detD' 2 ) 



-P.(2) 



(C4) 



where a.-\j is the residue of resonance i = 1,2 with pole at z — Zqj. 

Using Eqs. (O and ([5]), the pole residues r = \r\e' e as quoted by the PDG l6lll can be 
calculated. For the residue phase 9 16111 we consider the usual II 1 3L 1 1 4fl definition given by 



T = T B + 



M ■ 



■ iT/2 



(C5) 



for a resonance with width F on top of a background Tg. Comparing Eq. (lC.5t with Eq. <j3j and 
using Eq. (O, the pole residue r and its phase are given by 



\a-\PnNV 



-Ti + arctan 



Im (a-\p nN ) 



Re (a_i p nN ) 



(C.6) 



where p„n is the phase space factor p from Eq. (O for the nN — > nN transition, evaluated at the 
complex pole position. For the corresponding quantity in the reaction n + p — > K + T, + , one simply 
replaces p„ N -> y/p KN p K1 _. 

It is convenient to express the n 2 different residues [with i, f = 1, ■ ■ • , n for the transi- 
tions within n channels] in terms of n parameters g. Indeed, the residues factorize with respect 
to the channel space and, e.g., for the residues into the nN and KI, channels, 



a l -i =gigf 



(C7) 



with a unique set of g, up to one undetermined global sign. For the channels nN and Kl., the 
partial decay widths are evaluated using 



\gj\ 2 Mf l^fEf 
Y f (M R ,M f ,m f ) = k, g f = 2ng f 



In M R 



M 



(C8) 



/ 



with the resonance mass (final baryon, final meson mass) Mr = Rezo (Mf, mj) and u)f, Ef are 
the meson and nucleon energy at the on-shell momentum k. 
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The sum of partial decay widths should equal the total width, D/T/ = r tot . The right- 
hand side of this equation can be determined independently (T tot = -2/Imzo) and be used as 
a test of the formalism. Indeed, below the nnN threshold, the equality holds to the 1 % level 
(The definition of branching ratios into the effective nnN states pN, crN, nA is model-dependent 
anyways). Although Eq. ( IC.8I 1 is a good approximation to the partial decay widths, it should be 
noted that 1/ = r tot never holds exactly, even in a manifestly unitary coupled-channels model 
with only stable intermediate states. This is simply because the amplitude has non-analytic 
branch points, required by unitarity, and this information is not contained in the g. However, this 
does not become a real issue unless a pole is very close to a branch point. 
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